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We study the time-dependent Dalitz plot of D — > Ksiv + n~ in B° —> D^h decays, where h° is 
a 7T°, 77, rj , or u; meson and D* — » D/, using a data sample of 383 x 10 6 T(4S) — > BB decays 
collected with the BABAR detector. We determine cos 2/3 = 0.42 ± 0.49 ± 0.09 ± 0.13, sin 2/3 = 
0.29±0.34±0.03±0.05, and |A| = 1.01±0.08±0.02, where the first error is statistical, the second is 
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the experimental systematic uncertainty, and the third, where given, is the Dalitz model uncertainty. 
Assuming the world average value for sin 2/3 and |A| = 1, cos 2(3 > is preferred over cos 2/3 < at 
86% confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Time-dependent CP asymmetries in B° meson decays, 
resulting from the interference between decays with and 
without B°-B° mixing, have been studied with high pre- 
cision in b — > ccs decay modes by the BABAR and Belle 
collaborations These studies measure the asymme- 
try amplitude sin 2(3, where (3 = -axg(V c< iV* b /VtdV* b ) 
is a phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix [2], V . The CP violating phase 2/3, 
inferred from sin 2/3, has a two-fold ambiguity, 2/3 and 
■k — 2(3 (four- fold ambiguity in (3). This ambiguity can 
be resolved by studying decay modes that involve multi- 
body final states B° -► J/ipKhr° Q, D[K° s ir+ir~}h° [J, 
D*+D*~K° @ or K + K~K° Q, where the knowledge of 
the variation of the strong phase differences as a function 
of phase space allows one also to measure cos 2/3. 

In this Letter, we present a study of CP asymmetry 
in B° — > D^h° 0] decays with a time-dependent Dalitz 
plot analysis of D — ► K^tt + tt~ Q, where h° is a 7r°, 77, 
r]', or u> meson. The B° — > D^h° decay is dominated 
by a color-suppressed b — * cud tree amplitude. The di- 
agram b — > ug?, which involves a different weak phase, is 
suppressed by V u bV*JV c bV* d ~ 0.02. Neglecting the sup- 
pressed amplitude, we factorize the decay amplitude of 
the chain B° -> ZJ°/i° -> [i^Tr+Tr"]/! into .4/ = „4 B „4; 



^ g ^ D o. The D u 



and D' 



and similarly for B into A 

decay amplitudes are functions of the Dalitz plot vari- 
ables Ado — f{m 2 + ,m 2 _) and A^ = /("t-h 



/(m 2 



, to; 



-) 



where m\ 



._. In the T(45) -> 



system, the rate of a neutral 5 meson decaying at proper 
decay time i rec , the other B (Stag) at ttag> an d the D de- 
caying at a point on the Dalitz plot, is proportional to 



-FAt 



-\A L 



T{\A E0 \ 
±2|A|^o(-l) L Im(e 



(|^„| 2 + |A| 2 |^o| 2 ) 
- |A| 2 |^ D o| 2 )cos(ATOAi) 



(1) 



' 2 ^^o^ )sin(ATOAi) 



where the upper (lower) sign is for events with £?t ag de- 
caying as a B° (B°), At — t ICC —tt ag , T is the decay rate of 
the neutral B meson, A = e~ 2l ^{Ag/ Ab), Am is the B°- 
B° mixing frequency, £^0 is the CP eigenvalue of h , and 
(— 1) L is the orbital angular momentum factor. Here we 
have assumed CP-conservation in mixing and neglected 
decay width differences. For Dh° modes, £/i ( — 1) L = — 1. 
For D*[DTr°}h° (h° =1 uj) modes £ h o(-l) L = +1 including 
factors from D* decay [9j . In the last term of expression!]] 
we can rewrite 

lm(e' 2ll3 A D oA* no ) = lm(A D oA* Bo ) cos 2/3 

-Re(^o^ )sin2/3, (2) 



and treat cos 2(3 and sin 2/3 as independent parameters. 

We fully reconstruct B° — » D^h candidates from a 
data sample of (383±4) x 10 6 T(45) decays into BB pairs 
collected with the BABAR detector at the asymmetric- 
energy e + e _ PEP-II collider. The BABAR detector is 
described in detail elsewhere The decay modes 

used are Dir°, Drj, Drf , Duo, D*tt°, and D*-q, with 
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D* -► Dtt°, D -> K°ir + n-, K° s -► tt+tt", tt c 
jy — j. yy, 7r°7r + 7r _ , 77' — * ?7 7r + 7r _ , and — * 7r°7T + 7r~. 

Charged tracks are considered to be pions. The ifg 
candidate is reconstructed from 7r + 7r~ pairs, whose \ 2 
probability of forming a common vertex is greater than 
0.1%, with invariant mass within 10MeV/c 2 of the nom- 
inal K® mass The distance between the K® decay 
vertex and the primary interaction point projected on the 
x-y plane (perpendicular to the beam axis) is required to 
be greater than three times its measurement uncertainty. 
The angle 9k between the K® momentum and the line 
connecting the production and decay vertices of the K® 
on the x-y plane is required to satisfy cosOk > 0.992. 

An energy cluster in the electromagnetic calorimeter, 
isolated from any charged tracks and with the expected 
lateral shower shape for photons, is considered a photon 
candidate. A pair of photons forms a ir° — > 77 (77 — > 77) 
candidate if both photon energies exceed 30 (100) MeV 
and the invariant mass of the pair is between 100 and 
160 MeV/c 2 (508 and 588 MeV/c 2 ). If the 77 is paired with 
a D* , the invariant mass window is tightened to 515 < 
to 77 < 581 MeV/ c 2 . The 77 — » 77 candidate is rejected 
if either photon, when combined with any other photon 
in the event, has an invariant mass within 6 MeV/c 2 of 
the nominal tt° mass. We perform a kinematic fit to the 
photon pair with its invariant mass constrained at the 
nominal 7r° or 77 mass and reject candidates with a fit 
probability less than 0.1%. 

The t]/uj — > 7r°7r + 7r~, 77' — * rjit + it~, and D — > 
KgTT + ir~ candidates are formed by combining a ir°, 77, 
or K® with two charged pions. The % 2 probability of 
the decay products coming from a common vertex for 
h° (D) is required to be greater than 0.1% (1%). The 
momentum of the 7r° and the 77 candidates used in oj 
and 77' reconstruction must be greater than 200MeV/c. 
The invariant masses of the 77, 77', and lo candidates are 
required to be within 10, 8 and 18 MeV/c 2 of their respec- 
tive nominal masses, which correspond to approximately 
twice the RMS of the signal distributions. We retain D 
candidates within 60 MeV/ c 2 of the nominal D° mass, ap- 
proximately 10 times its mass resolution, to include suffi- 
cient data in the sideband. A kinematic fit is performed 
on the D candidate to constrain its mass to the nominal 
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D° mass. A D* — > -D-7T candidate is accepted if the in- 
variant mass difference between D* and D candidates is 
within 3 MeV/ c 2 of the nominal mass difference. 

The signal is characterized by the kinematic variables 
toes = V(s/2 + Po • Pb) 2 /^o - Pb» and A£ = £* - 
^beam' wnere the asterisk denotes the quantities evalu- 
ated in the center-of-mass (cm.) frame, the subscripts 0, 
beam and B denote the e + e~ system, the beam and the 
B candidate, respectively, and yfs is the cm. energy. For 
signal events, toes peaks near the B° mass with a reso- 
lution of about 3MeV/c 2 , and AE peaks near zero, with 
a resolution that varies by mode. We require uies > 
5.23GeV/c 2 and select events with \AE\ < 80 MeV for 
modes with 7r°, 77 — > 77, and \AE\ < 40 MeV for modes 
with rj,uj — ► 7r°7r + 7r _ , or rf — * ^7r + 7r _ . 

The proper decay time difference At is determined 
from the measured distance between the two B decay ver- 
tices projected onto the boost axis and the boost (/J7 = 
0.56) of the cm. system. The reconstructed |Ai| and its 
uncertainty a At are required to satisfy |At| < 15 ps and 
a At < 2.5 ps. The flavor of Stag is identified from parti- 
cles that do not belong to the reconstructed B meson us- 
ing a neural network based flavor-tagging algorithm [l2T ] . 

The main background is from the continuum e + e~ — > 
qq (q = u,d,s,c). We use a Fisher discriminant (J 7 ) 
to separate the more isotropic BB events from more 
jet-like qq events [1 31 ) - The requirement on T is op- 
timized with simulation. Another major background 
for the D*tt° mode comes from color-allowed B~ — > 
D°p~(p~ — > 7r°7T~) decays, which mimic signal if the 
7r~ is missed from reconstruction while a random tt° is 
included. We veto the B° candidate if the combination of 
another charged pion in the event with the D and the ir° 
in the B° candidate is consistent with a charged B decay. 
In total we select 4450 events, of which 2843 events have 
useful tagging information (tagged). 

The signal and background yields are determined 
by a fit to the (mEs,wiD) distributions using a two- 
dimensional probability density function (PDF), where 
rriD denotes the K®n + ir~ invariant mass. We divide the 
sample into four categories to take into account different 
background levels: (1) Dir°, (2) Drj and Drf (3) Doj, 
and (4) D*h°. The PDF has five components: (a) signal, 
and backgrounds that peak in (b) both Toes an d to 73 , ( c ) 
Toes but not mo, (d) mjj but not toes, and (e) neither 
distribution. Both peaks are modeled by a Crystal Ball 
line shape [141 ]. The non-peaking component is modeled 
by a straight line in me and a threshold function [15} 
in Toes- We fit the four categories of events simultane- 
ously, allowing the toes peak shape to be different but 
letting them share the tud shape and toes background 
parameters. We first determine the amount of the peak- 
ing component (b) from simulated events and then fit to 
data allowing all other components to vary. We obtain 
463 ± 39 signal events (335 ± 32 tagged) . The contribu- 
tion from each mode is shown in Table |U The Toes and 
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FIG. 1: Distributions of (a) m E s [\m D -m P I % G \ < 14MeV/c 2 ], 
and (b) D mass [tties > 5.27GeV/c 2 ]. Dashed (dotted) lines 
represent the total (non-peaking) background. 



mo distributions are shown in Fig.[TJ 

The D° — > K®ir + ir~ Dalitz plot has been studied in 
detail 

HE E3]. We 

use the isobar formalism described 
in [l8| to express Apo as a sum of two-body decay matrix 
elements {A r ) and a non-resonant (NR) contribution, 



A d q = a NR e 



A r (m 2 + ,m 2 _). (3) 



The function Arim 2 ^, m 2 _) is the Lorentz-invariant ex- 
pression for the matrix element of a D° decaying into 
K®7r + Tr~ through an intermediate resonance r, parame- 
terized as a function of the position on the Dalitz plot. 
The resonances in the model include if* (892), K * (1430), 
#1(1430), AT* (1410) and if* (1680) for both K° s tt+ and 
K° s n-, and p(770), c^(782), / (980), / (1370), / 2 (1270), 
p(1450), and two scalar terms a and a' in the ir + ir~ sys- 
tem. Details of the Dalitz model and the parameters 
(determined from data) can be found in [rH |. 

To perform the time-dependent Dalitz plot analysis, 
we expand the PDF to include At and Dalitz plot de- 
pendence. The signal component is proportional to ex- 
pression [TJ modified to account for the probability of 
mis-identifying the Stag flavor (mistag), and is convolved 
with a sum of three Gaussian distributions 19]. The 
mistag parameters and the resolution function are deter- 
mined from a large data control sample of B° — * D^*^h + 
decays, where h + is a 7r + , p + , or meson. Each of 
the background components consists of a product of At 
and (m+jiriL) PDFs. The components that peak in 
use Aijo{m 2 ,^rn 2 _) as their Dalitz model. The model for 
components that are flat in mo is an incoherent sum of 
a phase space contribution and several resonances. The 
choice of resonances and their relative contributions are 
determined empirically from events outside the mo peak. 
The At model for components that peak in tties is a 
simple exponential decay convolved with the resolution 
function used in the signal component. For the non- 
peaking background, we use a zero-lifetime component 
convolved with a double- Gaussian resolution function for 
events with a real D because they are dominated by cc 
events, and we add an exponential decay component for 
events without a real D to account for B background. 
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FIG. 2: Dalitz plot projections for (a,b,c) _B°-tagged events 
and (d,e,f) _B°-tagged events. Points are data, open his- 
tograms are PDF projections, and shaded histograms are 
background contributions. 



We fit the mES> fnr>, and At distributions, with mES 
and mo shapes and background fractions fixed by the 
previous fit for event yields, to determine the At pa- 
rameters for backgrounds. We then perform the final fit 
adding Dalitz plot variables to determine cos 2/3, sin 2/3 
and |A|. Table Q] shows the nominal fit result (All) and 
the results of a fit allowing cos 2/3 and sin 2/3 to be dif- 
ferent among the four types of events. The correla- 
tions are p(cos 2/3, sin 2/3) = 2%, p(|A|, cos 2/3) = 2%, 
and p(| Aj, sin 2/3) = —2%. The Dalitz plot projections 
are shown in Fig. [5] Figure [3] shows the time-dependent 
asymmetries (N + - N-)/(N+ +iV_), where N+(N-) is 
the number of B°(B°) tagged events, for events in vari- 
ous Dalitz plot regions. Events in the D — > K®p region 
are dominated by a single CP eigenstate, thus the asym- 
metry is proportional to sin 2/3 sin(AmAt). Events near 
D — > i4T* ± 7r =F are dominated by decays to a definite fla- 
vor, and therefore exhibit a cos(AmAt) behavior. 



TABLE I: Tagged event yields iV tl 
statistical. 



and fit results. Errors are 



Mode 


AW 


cos 2/3 


sin 2/3 


|A| 




143 ±19 


0.78 ±0.92 


0.70 ±0.52 




DnM 

Duj 


60±f2 
76 ±12 


1.20 ±1.19 
0.43 ±0.87 


-1.17±1.00 
-0.48 ±0.74 


1.0 (fixed) 


D*h° 


56 ±12 


-0.56 ±1.07 


0.78 ±0.87 




All 


335 ±32 


0.42 ±0.49 


0.29 ±0.34 


1.01 ±0.08 



The dominant systematic uncertainty is the Dalitz plot 
model dependence. The Dalitz model includes scalar 
terms a and a' which are not well established, in or- 



der to achieve a good quality fit [17| • We study the effect 
of these two scalars by simulating a number of datasets, 
each of which is 50 times the size of the data, accord- 
ing to the PDF, and repeat the final fit using both the 
nominal PDF and the PDF without the two scalars. We 




5 -5 
At (ps) 



5 
At (ps) 



FIG. 3: Time-dependent asymmetries for (a) D — > K%p region 
(\m n + n - — 770| < 150), where the opposite CP asymmetry in 
D*h° has been taken into account, (b) D — > K* + ir~ region, 
and (c) D — » K*~ix + region (\m K o^ — 892 1 < 50). Units are 
MeV/c 2 . Curves are projections of the PDF. 



compare the results between the two fits in each dataset 
and conservatively take the quadratic sum of the mean 
and RMS of the differences as the systematic uncertainty: 
cr(cos2/3) = 0.13, cr(sin2/3) = 0.05, and a(\X\) < 0.01. 
Many parameters are pre-determined in fits to control 
samples and to data without the Dalitz variables. We 
randomize them according to a Gaussian distribution 
whose width equals one standard deviation of each pa- 
rameter, taking correlations into account, and repeat the 
final fit. The width of the distribution is taken as the 
systematic uncertainty: er(cos2/3) = 0.06, <r(sin2/3) = 
0.02 from Dalitz model parameters; c(cos2/3) = 0.05, 
cr(sin2/3) = 0.02 from and wes shape parameters; 
cr(cos2/3) < 0.01, cr(sin2/3) < 0.01 from background At 
parameters, tagging parameters, or signal At resolution 
function. We also vary the peaking background fractions 
by the statistical uncertainty found in simulation and find 
the variations are cr(cos2/3) = 0.02 and cr(sin2/3) = 0.01. 
Other sources of uncertainty such as B°-B° mixing fre- 
quency, B lifetimes, background Dalitz model and re- 
construction efficiency variation over the Dalitz plot are 
negligible. In all cases, the uncertainty on |A| is less than 
0.01. The only significant uncertainty on |A| (~ 0.02) 
is from the interference between the CKM-suppressed 
b — > ucd and CKM-favored b — > cud amplitudes in some 
Stag final states [20(. This effect is studied with simu- 
lation. Summing over all contributions in quadrature, 
we obtain total experimental systematic uncertainties 
cr(cos2/3) = 0.09, cr(sin2/3) = 0.03, and ct(|AQ = 0.02. 

To resolve the ambiguity in 2/3, we generate two sets 
of toy simulation samples, one with cos 2/3 = \J\ — S\ = 
Co, and the other with cos 20 = —Co, where So = 0.678, 
the world average of sin 2/3 [2l| , and fit each sample while 
fixing sin 2/3 = So and |A| = 1. For data, this config- 
uration results in cos 2/3 = 0.43 ± 0.47. We then use 
double- Gaussian functions, h±{x) for ±Co hypotheses, 
to model the probability density of the resulting cos 2/3 
distributions, smeared by the experimental systematic 
uncertainty and the uncertainty of Co- The confidence 
level (C.L.) of preferring cos 20 = +Co over —Co is de- 
fined as h + (x)/[h + (x) ± h-(x)] if cos 2/3 = x is observed 
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in data. Considering the Dalitz model dependence for 
cos 2/3 (0.13), we use x between 0.43 ± 0.13 and find the 
smallest C.L.= 86% at x = 0.43 - 0.13. 

In conclusion, we have studied the B° -> D^h° 
decays using a timc-dcpcndcnt Dalitz plot analysis of 
D -> K^tt+tt-. We obtain cos 2/3 = 0.42 ± 0.49(stat.) ± 
0.09(syst.) ± 0.13(Dalitz), sin 2/3 = 0.29 ± 0.34(stat.) ± 
0.03(syst.) ± 0.05(Dalitz), and |A| = 1.01 ± 0.08(stat.) ± 
0.02(syst.). Using the world average sin 2/3 = 0.678 ± 
0.026 and |A| = 1, cos2/3 > is preferred over cos2/3 < 
at 86% C.L. 
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